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Digital video-microscopy measurements are reported of both elastic bandstructures and over- 
damped phonon decay times in two-dimensional colloidal crystals. Both quantities together allow 
to determine the friction coefficients along various high symmetry directions in <f-space. These 
coefficients contain valuable information on the hydrodynamic forces acting between the colloidal 
particles. We find Stokes- like friction for phonons near the edge of the first Brillouin zone and van- 
ishing friction coefficients for long wavelength phonons. The effect of this wavelength dependence in 
real-space is further investigated by simulating a crystal with constant friction (Langevin simulation) 
and comparing experimentally measured and simulated particle auto-correlation functions. 



I. INTRODUCTION 

The lattice dynamics of crystals composed of colloidal 
particles is entirely different from that of atoms. Being 
immersed in a viscous solvent, colloids experience fric- 
tion which strongly dampens their motion. As a result, 
phonons in colloidal systems - characterized by the po- 
larization index j and the wave- vector q - show an over- 
damped dynamics and decay exponentially with decay 
times T(qj). These times are related to the eigenvalues 
X(qj) of the dynamical matrix which are refered to as the 
elastic bandstructure in the following. These eigenval- 
ues may be pictured as the "spring-constant" associated 
with the phonon (qj). In overdamped systems, the rela- 
tion between T(qj) and X(qj) is the phonon dispersion 
curve which replaces the relation between the frequency 
of the propagating phonons and the bandstructure X(qj) 
in atomic systems. The phonon dispersion curve quanti- 
fies the response of a crystal to perturbations of a given 
wavelength; as such it is the central relation of every crys- 
tal, characterizing both its static and dynamical proper- 
ties. 

Contrary to atomic systems, where the phonon disper- 
sion curve is entirely determined by interparticle forces, 
in colloidal systems hydrodynamic interactions have to 
be considered. These interactions arise when moving col- 
loids exchange momentum through the solvent. Because 
of their long-ranged nature these interactions are difficult 
to treat theoretically, especially in crystals pQ. It has 
been demonstrated by Hurd et al. [2j [3] that along cer- 
tain high symmetry directions within the colloidal crys- 
tal hydrodynamic interactions can be taken into account 
through the wavelength dependence of the friction co- 
efficients "f(qj). These coefficients connect the phonon 
decay times and the bandstructure such that in the limit 
of strong damping ("j(qj) S> X(qj)) the phonon disper- 
sion relation becomes T(qj) = "f(qj)/ X(qj). Hence, with 
the wavelength dependence of the friction coefficients one 
obtains detailed information also on the hydrodynamic 



forces acting within colloidal crystals. 

The purpose of this paper is to report measurements 
of colloidal crystal dispersion curves of thermally excited 
overdamped phonons, and, more specifically, measure- 
ments of the (f-dependence of the friction coefficients in- 
volved. We will compare the experimental data with 
Langevin simulations where the friction is taken into ac- 
count only by the constant Stokes friction. This com- 
parison will demonstrate what effect the wavelength de- 
pendence of the friction can have on correlations in real 
space. The simulations also serve to check for the consis- 
tency of the evaluation procedure, for finite size effects 
as well as possible sampling errors. 

Measuring cf-dependent friction coefficients in colloidal 
crystals is experimentally demanding, and such experi- 
ments have - to our knowledge - never been attempted. 
Such studies require the simultaneous measurement of 
both the decay times and the bandstructure. While decay 
times are experimentally accessible with dynamic light 
scattering techniques [51 |U [5J [7] or Brillouin scattering 
[HI H], the experimental determination of X(qj) requires 
real space positional information. Due to the mesoscopic 
time- and length scales, such information is accessible in 
colloidal systems with digital video microscopy which al- 
lows to follow particle positions with a resolution of a few 
nm at video rates as high as 25 Hz. 

The literature on phonon dynamics in overdamped sys- 
tems is still manageable. Of central importance to the 
present paper is the work of Hurd et al. [2] who starting 
from [1] not only present a theory of hydrodynamic in- 
teraction in colloidal crystals but also performed photon- 
correlation spectroscopy measurements to study the dis- 
persion of lattice waves. Early studies of single-particle 
dynamics within colloidal crystals employed scattering 
techniques [TDJ[TT], while the first video-microscopy ex- 
periments relevant for our question were presented by 
Bongers and Versmold [12] who measured the particle 
mean square displacement as a function of time, later the- 
oretically analysed within the framework of the Langevin 
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model in (15) [II]. 

Hydrodynamics within colloidal systems that are not 
in the crystalline state, is the subject of a large number 
of papers. Digital video microscopy has been used to in- 
vestigate dynamic properties (dynamic structure factor, 
the hydrodynamic function, hydrodynamic diffusion co- 
efficients) in colloidal suspensions confined between two 
parallel glass plates [T5J [TBI [T7j. Colloidal friction has 
also been a major issue in recent optical tweezer mea- 
surements of the hydrodynamic interaction between two 
colloidal particles [TH1 EH HD] ■ The time-independent hy- 
drodynamic forces as well as the cross-correlation func- 
tion of the colloid position have been interpreted in terms 
of two-body hydrodynamic interactions, embodied by the 
Oseen tensor. However, a two-body description of hydro- 
dynamics in crystals is inadequat because of the neglect 
of many-body effects which dominate the hydrodynamic 
interaction in regular arrays of particles [TJ [5] . 

When at long wavelengths the frictional force vanishes, 
damping can become so weak that even in colloidal sys- 
tems phonons should be able to propagate. Hurd et al. 
predicted a switch to such propagating modes for trans- 
verse modes but were not able to observe them, proba- 
bly because of the finite-cell geometry and hydrodynamic 
wall effects - a conclusion that has later been corrobo- 
rated by Derksen et al. |4J who carefully reconsidered the 
Hurd experiment. An alternative explanation offered by 
Felderhof and Jones [2T] in terms of additional damping 
through retarded counterions was challenged by Derk- 
sen et al. [1] and Hoppenbrouwers et al. [5]. Evidence 
for overdamped transverse modes turning propagative at 
long wavelengths has later been produced by Tata et al. 
7 in finite size millimeter crystals. We here completely 
ignore such propagating modes, returning to a brief dis- 
cussion of their possible existence only in the last section. 

The paper starts with a short presentation of the over- 
damped Langevin model of a colloidal crystal, which 
forms the basis of our simulations. Then, the experi- 
mental setup is described, followed by a section in which 
we outline how bandstructures, friction coefficients and 
correlation functions are obtained from particle configu- 
rations. We proceed with the coefficients and the parti- 
cle auto-correlation functions to show the differences be- 
tween systems with and without wavelength-dependent 
friction, i.e., the differences between simulation and ex- 
periment. The main conclusions of the paper are formu- 
lated in the last section. 



II. EXPERIMENTAL 

A. Simulation: The overdamped Langevin model 
of a colloidal crystal 

We consider a two-dimensional (2D) hexagonal crystal 
of N colloidal spheres with diameter a immersed in an 
aqueous electrolyte being characterized by the dielectric 
constant e and the inverse screening length n. Particle 



positions at time t are denoted by x n (t) with n label- 
ing a hexagonal lattice site R n . The position vector de- 
composes into x n (t) = R n + u n {t) with u n {t) being the 
particle displacement from the lattice site R n . The inter- 
action of two colloids at distance r is given by the Yukawa 
potential, 



$(r) 



(Z e ffe) 2 exp(Kcr) exp(— nr) 
Ane e (1 + kct/2) 2 r 



(1) 



where Z e ff is the effective colloidal charge. Within the 
framework of the harmonic approximation, the particles 
interact through elastic forces characterized by the spring 
constant /cq, given by 
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where a is the nearest-neighbor distance in the hexagonal 
crystal. 

As explained in the next section, in our experiments 
the colloidal crystal is stabilized by a commensurate 
hexagonal, light-induced substrate created by three in- 
terfering laser beams. In the sample plane, the electric 
fields of the linearly polarized laser beams are given by 

4(f) = l exp(i^x) i = 1,2,3 (3) 
Kae,, , 



with K\ 
K 2 
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where e x — [1, 0] and e y — [0, 1] are the basis vectors 
of Cartesian coordinates, and K is chosen to match com- 
mensurate conditions. The interfering laser beams create 
the substrate potential 



^2 exp(iKix) 
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4H3-^cos [Kijx] ) 
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where all three vectors Kij — Ki — Kj must have a length 
K' = 47t/(v / 3a) to ensure a fully commensurate sub- 
strate. To adapt this potential to the harmonic approxi- 
mation used in our simulations, we expanded the cosine 
function in eq. Q for positions x n = R n + u n near a 
lattice site R n , and finally obtain 



$ext{Rn + U n ) 



A 2 a K' 2 
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where = K^jK' and where use was made of Kij ■ 
R n = Iran with an integer m. Thus, we arrive at another 
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spring constant 



B. 



fci := A 2 Q K' 2 = 



A 4n 
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(0) 



denning the strength of the springs with which every par- 
ticle is pinned to its lattice site R n . 

Having introduced k and fcj we can now proceed to 
the model equation describing the dynamics of the col- 
loidal crystal. We here adopt a simplified Langevin de- 
scription of the lattice dynamics, ignore hydrodynamic 
interactions and include just the Stokes friction between 
the colloids and the solvent. For more elaborate lattice 
dynamic theories for colloidal crystals including hydro- 
dynamic interactions, see |H1 11H HI HI • 

If R n is an arbitrary reference site of the hexagonal 
lattice with its six nearest neighbors R\ = R n +ae y , i?2 = 
Rn — \[^>aj2e x + a/2e y , and so forth, the overdamped 
Langevin equation in the harmonic approximation reads 
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du na (t) 



dt 
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^ ^ (Rma — Rna){Rmf3 — Rnfi) 
X [u m (t) - U n (t)} + kiU na (t) + f na (t) = 0. 



(7) 



where we have restricted ourselves to nearest-neighbor in- 
teractions and where an insignificant term ~ 5 a p<fr'(r)/r 
has been omitted. f n (t) is a randomly fluctuating force 
for which the following thermal averages must be satisfied 



{fna(t)fn 



(fna(t))=0, 

,(t + r)) =2 1 k B T5(T)S n 



i' $ac 



(8) 



with fcrjT being the thermal energy. In eq. Q, we have 
introduced the Stokes friction coefficient 7 which mea- 
sures how strongly a quiescent fluid resists the motion of 
an isolated sphere. Clearly, this is a rough approxima- 
tion as the results of the present paper will once more 
confirm: neither can a crystal particle be considered iso- 
lated, nor is the surrounding fluid quiescent, but will in 
fact be itself in a dynamical state due to the motion of 
all other particles. 

One can express the times in units of 7 pm 2 /k B T and 
the spring constant kg in terms of feeT/pm 2 such that 
the number of independent input parameters to eq. ^ is 
reduced to just three values. In our simulations we have 
chosen k \±m 2 /k B T — 20, fei/fco = 0-6 and a = 4 pm. 
While the simulations have been performed with reduced 
units, the data will be presented in normal physical units. 

Finally, it must be mentioned that for technical reasons 
a value of k [im 2 /k B T = 20 and above turned out to be 
the most convenient choice. These values are well above 
the ones observed in our experiments; however, the band 
structure scales with fco as has been shown in [25] and, 
therefore, our simulations allow for a meaningful analysis 
of the experimental results. 



Experiment: 2D colloidal crystals stabilized by 
light-induced substrates 




FIG. 1: Sketch of the sample cell containing the 2D crystal 
subject to a light-induced substrate (indicated for one parti- 
cle) and, additionally, pressed down towards the bottom plate 
by an normal incident laser beam. 

Experiments were performed with an aqueous suspen- 
sion of highly charged polystyrene spheres with a = 
2.4 pm diameter and a polydispersity below 4%. The par- 
ticles interact via a screened Coulomb potential, eq. (JlJ , 
with a renormalized surface charge of about Z e ff « 
14000 [26 and a screening length ss 458 nm, de- 
termined as described further below. As sample cell we 
used a cuvette made of fused silica with 200 urn spacing 
between top and bottom plate which was connected to 
a standard closed deionization circuit to maintain stable 
ionic conditions during the measurements. An Argon-ion 
laser beam (wave length = 488 nm) was scanned around 
the central region of the sample to create a boundary 
box whose size could be continuously adjusted by a pair 
of computer-controlled galvanostatically driven mirrors. 
This allowed us to adjust precisely the particle area den- 
sity [27]. At sufficiently high particle densities the parti- 
cles form a hexagonal 2D colloidal crystal (a = 6 pm) 
close to the bottom plate (see Fig. [IJ. Out-of-plane 
fluctuations are reduced by pushing the particles down 
towards the bottom plate with a normal incident laser 
beam. Due to the interplay of the electrostatic particle- 
wall repulsion and the light pressure from above, this re- 
sults in a mean particle-wall distance of about 2 pm and 
normal fluctuations of an amplitude of less than 100 nm. 
By systematic variation of the light pressure we have ver- 
ified that within our parameters it does not influence the 
results discussed below. 

It is known that 2D crystals cannot have perfect long- 
ranged translational order because low-<f phonons desta- 
bilize the system. This would considerably reduce the 
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quality of our measured particle configurations and the 
correlation functions generated from it. We therefore sta- 
bilized the crystal by exposing it to the perfectly hexago- 
nal commensurate substrate potential from eq. Q being 
created from the interference of the three laser beams 
in eq. § (P = 5W, frequency-doubled Nd:YV04 laser, 
wavelength = 532 nm). It has been shown theoretically 
and experimentally J25[ US] that if the light-induced sub- 
strate is created from laser beams of the same intensity 
the resulting bandstructure is identical to that of the free 
crystal, except that it is shifted by ki. 

Particle positions and trajectories were determined for 
a 16 shell hexagonal crystallite (N = 817 particles) which 
was part of a much larger crystal comprising about 2000 
to 3000 particles. All derived quantities were computed 
from averages over 29500 images, using digital video mi- 
croscopy at an aquisition rate of 8 frames per second. 



that along certain high symmetry directions (path sec- 
tions 1 and 3, see central graph in top row in Fig.[2|, hy- 
drodynamic interactions lead to an overdamped phonon 
dispersion relation of the more general form, 



T(qj) 



7(gj) 



(12) 



where the friction coefficient ^(qj) is no longer constant, 
but depends on the band index and the wave vector q. 

Having introduced these equations, we can formulate 
more precisely the central idea of the present work. Equa- 
tion (pi allows us to determine the bandstructure by tak- 
ing X[qj) = k-e,T/A{qj). Equation (10 1, on the other 



hand, permits an independent determination of the decay 
times T(qj). Multiplication of both quantities according 
to eq. ( 12 1 gives us access to the wave- vector dependent 



friction coefficients, 



C. Data processing: from configurations to derived 
quantities 

Simulation and experiment produce the same output: 
sets of particle configurations. These will then be fur- 
ther processed to generate what may be called "derived 
quantities" . Both sets of configurations, the computer- 
generated ones as well as the experimental ones, will be 
processed in the same way. 

It is a straightforward procedure to compute for a 
given configuration the normal coordinates Q{qj, t) di- 
rectly from the measured displacements u n (t) [25]. Here, 
q is one of the N allowed phonon wave vectors, while j 
is the branch index, j = 1,2, distinguishing between the 
longitudinal and transversal branch. The normal coordi- 
nates allow us to determine the normal mode bandstruc- 
ture X(qj) through the relation 



Am = (Q*mQ(m = 



k B T 



(9) 



Several experimental bandstructures of colloidal systems 
being based on this relation can be found in literature 
[28} |29"1 EOj. Starting from the overdamped Langevin 
equation d7|, one finds that the normal-mode auto- 
correlation function decays exponentially 



(Q*(qj, t + r)Q(qj, t)) = A(qj) exp - 



T(qj)) 



with the overdamped phonon dispersion relation 
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T(qj) = 



Hqj) 



(10) 



(ii) 



connecting the phonon decay times T(qj) with the elastic 
bandstructure X{qj). So far, the friction coefficient 7 has 
been taken to be constant, relying on the assumption of 
the Langevin model. Going beyond this model, this latter 
equation takes a different form. Hurd et al. 2 showed 



i(qj) = X(qj)T(qj) 



(13) 



allowing a detailed experimental study of 7(<jj) and thus 
also a check of the Langevin model and its assumption 
of constant friction. 

The colloidal crystal dynamics can also be studied by 
examining the behavior of the particle auto-correlation 
function 



c(r) = ~ ^2(u na (t + T)u na {t)) 



N 
N ^ 



(Q(qj,t + T)Q*(qj,t)) 



(14) 



which is connected to the time-dependent mean-square 
displacement 5msd(t) through 

C ( T ) = ^msd(oo) - <5 M sd(t)) ■ (15) 

Within the Langevin description, this auto-correlation 
function is related to the Laplace transform of the phonon 
spectrum [25] 



fcoc(r') 



dX'^p-exp{-\'T') (16) 



k B T J X 
where r' = rfco/7, A' = X/k and where 



(17) 



is the phonon spectrum, counting the number of X(qj) 
within the interval [A, X + dX]. This relation takes a static 
quantity - the system X(qj) of (^-dependent spring con- 
stants - to predict a quantity c(t) that characterizes the 
dynamics of the system. In classical atomic crystals, a 
similar relation exists, connecting the Fourier transform 
of the velocity auto-correlation function with the phonon 
density of states [31] . Equation (16 1 is, of course, valid 



only within the framework of the Langevin model. 
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III. RESULTS 

A. bandstructure, phonon decay and friction 
factors 

Figure [2]shows phonon auto-correlation functions com- 
puted either from simulation data (left column of plots) 
or from experimental configurations (right column of 
plots) for four different allowed wave vectors q in the 
first Brillouin zone (center column of plots). As outlined 
in eqs. (10 1 to (13 1, we expect an exponential decay and 
therefore fitted the data to the function Aexp(— t/T). 
While such an exponential decay is indeed observed in 
the Langevin simulation data the experimental data show 
some deviations at larger times. The experimental and 
simulation curves are obtained from averages over 29500 
and 10000 configurations, respectively. Averaging over 
time windows of different lengths, one can show that for 
correlation functions near to the edge of the first Brillouin 
zone (BZ) these deviations of the experimental curves 
are due to insufficient sampling. However, for correla- 
tion functions near the center of the first BZ, one finds 
deviations which seem to be independent of the average 
procedure. This might suggest that these (rather small) 
deviations are fingerprints of propagating modes which 
are not considered in our overdamped Langevin descrip- 
tion and can thus not be reproduced with exponential fit 
functions. Such propagating modes are expected when 
the system is damped, but not overdamped, i.e., when 
the damping is not strong enough to prevent an oscilla- 
tory behavior. The alternative explanation - that these 
deviations result from the typical 2D instability against 
low q phonons - can be safely ruled out because the sys- 
tem is pinned to a hexagonal lattice by the light-induced 
substrate. 

Performing such fits for all allowed wave vectors along 
the irreducible path 1 to 3 (see top row in Fig. [2]), one 
obtains the function T(cfj), that is, the phonon decay 
times for both branches as a function of the wave vec- 
tor. These times are plotted in the central row of Fig. [3] 
for both the experimental as well as the simulation data. 
On the other hand, one can use equation Q to deter- 
mine the elastic bandstructure, plotted for both sets of 
data in the top row of Fig. [3] Harmonic lattice the- 
ory of 2D colloids in light-induced potentials [25] pro- 
vides us with explicit expressions for the bandstructure; 
they are fitted in Fig. [3] to the measured structure us- 



path ABC 



D 



longitudinal 3 
longitudinal 1 
tranversal 3 
tranversal 1 



7.17 -1.45 0.07 
6.72 -1.24 0.07 
1.7 1.16 0.28 -0.02 
1.7 1.28 0.29 -0.09 



ing kg and k\ as fitting parameters. For the band- 
structures produced from the simulated data, we obtain 
ki/k Q — 0.49 and k Q — 18.9 k B T/[im 2 in good agreement 
with the values actually used in the simulation, while 
from the experimental structures we derive fci/fco = 0.61 
and fco = 3.3 fceT/u-m 2 . With that we can estimate with 
eq. (j2| (asumming Z e ff = 14000) a value for 1/k of 
458 nm which is actually quite a reasonable value. As 
expected, due to the presence of the stabilizing substrate 
the bandstructures are slightly shifted upwards, with a 
shift that is directly proportional to the constant k\. 



Following eq. ( 13 1 we multiply the two quantities X(qj) 



TABLE I: Fit-parameters of the function f(x) — A + Bx + 
Cx 2 + Dx 3 used in Fig. [Z] 



and T(qj) to obtain the friction coefficients 7(cfj), plotted 
in the bottom row of Fig. [3] In the simulations we used 
the value 7 = 8.33 feeTs/urn 2 and a similar value 7 = 
8.85 fee Ts/|J.m 2 , independent of the wave- vector, comes 
out from the evaluation procedure. This demonstrates 
the internal consistency of our scheme and the correct 
implementation of our data evaluation procedures. More 
importantly, it shows (i) that our statistics (which in the 
experiment is even better) is sufficient, and (ii) that finite 
size effects are negligible. While in the experiment we 
analysed a 16 shell hexagonal crystallite being part of 
a much larger crystal, the simulation has been carried 
out with 1000 particles and analysed within a 12 shell 
hexagonal crystallite comprising 469 particles. Finite size 
effects would show up in strong variations of ~f(qj) near 
the center of the BZ - something we do not observe in 
Fig. [3] Since the experimental system size is considerably 
larger than in the simulation, we can be rather confident 
that also the experimental data are not polluted by finite 
size effects. 

With that we turn to our main result - the function 
7(<fj) derived from the experimental data (Fig. [3^). We 
observe pronounced variations of the friction coefficient, 
having roughly constant values only in one band and only 
very close to the edge of the BZ. Most importantly, for 
both bands 7(<fj) tends to zero for q — > 0. Almost the 
same function ~f(qj) is obtained when the renormalized 
simulated bandstructure is multiplied with the experi- 
mental decay times. This finding demonstrates that it 
is essentially the differences in T(qj) which lead to the 
differences in j(qj) between experiment and simulation. 
The different q — > behavior of the friction coefficients 
can then be traced back to the decay times in this limit 
being considerably longer in the simulation than in re- 
ality, implying that experimentally the long wavelength 
phonons are decaying much faster than the Langevin sim- 
ulation with its constant friction assumption suggests. 
This is an effect due to hydrodynamic interactions that 
shall be further discussed in the last section. 

Figure [4] shows again the measured friction factors: 
now the coefficients for both q directions, path 1 and 3 
(see top row in Fig. [2j, appear over the same axis. Also, 
the coefficients have been averaged over three equivalent 
directions in q^space. Note that consideration of all other 
directions is obsolete due to symmetry reasons. To facil- 
itate easy reproduction of our results we have fitted the 
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FIG. 2: Phonon auto-correlation function for the two band indices j = 1 (closed symbols) and j = 2 (open symbols) at certain 
q values in the first Brilloun zone, as determined from Langevin simulation data (A-D) and from video-microscopy data (E-H). 
The center column shows the 817 allowed phonon wave vectors q in the first Brillouin zone; the q value actually considered 
in the respective row is represented as a black open symbol, g-vectors along the irreducible path are printed as solid black 
symbols, gray color is used otherwise. The numbers in the top graph label the sections of the path around the irreducible part 
of the first Brillouin zone. The solid lines are exponential fits to the data with the resulting phonon decay times given in the 
legend of each figure. The curves for j = 2 (open symbols) are vertically shifted for clarity. 
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Simulation Experiment 
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FIG. 4: Measured friction coefficients ~y(qj) as shown in 
Fig. [3^ together with fitted polynomials (lines). Fit param- 
eters are listed in table [T] Upper two curves show the coeffi- 
cients for the longitudinal modes, while the lower two curves 
correspond to the transversal modes. Open symbols for a di- 
rection in <f-space along path 1 (see top row in Fig.|5|, closed 
symbols for direction along path 3. 



FIG. 3: Elastic bandstructure A(gj') (top pair of graphs), 
phonon decay times T(qj) (central pair of graphs), and 
wavelength-dependent friction coefficient 7(gj) (bottom row), 
plotted along the irreducible path I to 3 of the first Brillouin 
zone (see top row in Fig. [2J. Solid lines in the upper graph 
represent bandstructures calculated within the framework of 
harmonic lattice theory. 



curves to a polynomial, with the fitting parameters given 
in table |TJ We observe for q — > that while the fitted 
functions of the longitudinal branches go to zero, those 
of the two transversal modes extrapolate our results to 
a value slightly larger than zero. We expect that this 
extrapolation is not very meaningful as the real curve 
between q = and our first measured data point might 
be completely different because of the emergence of prop- 
agating phonons in this q regime which would destroy the 
basis of our determination of T(qj). 



We finally observe that the theoretical Stokes friction 
coefficient Qnr]{a/2) for our colloidal spheres in water is 5 
kftTs / \±m 2 , showing that our coefficients are of a reason- 
able magnitude (see Fig. [3] and Fig. |3J. This statement 
is supported by the observation presented in [32 where 
7 increases to 8 — 10 feeTs/iim 2 if the colloidal sphere is 
located approximately 2 urn above a wall. 



B. Phonon spectrum and mean-square 
displacement 

Instead of considering the dynamics of all 2N phonons 
individually, one may also study the colloidal dynamics 
more globally by examining the average over all phonon 
auto-correlation functions (Q(qj, t+r)Q* {qj 1 1)), thus ar- 
riving at the particle auto-correlation function c(r) de- 
fined in eq. (14 1. As an average over all phonons, it is 



perhaps more suited to serve as a general measure of the 
colloidal dynamics; physically it is directly linked to the 
time dependent mean-square displacement, eq. ((151. 



To see how the observed differences between the actual 
colloid dynamics and the Langevin dynamics manifest it- 
self with respect to the particle auto-correlation function 
c(t), a Langevin prediction of c(r) for the experimental 
system is required. Here, we can exploit the relation ( 16 1 



which assuming a Langevin system derives a c(r) directly 
from the phonon spectrum G(A) and thus essentially from 
the bandstructure X(qj). A c(r) thus predicted can then 
be confronted with the actual particle auto-correlation 
function obtained from averaging (u na (t + T)u na (t)) as 
given in eq. (14). This is the idea of Fig. [5] 

The figure shows almost identical phonon spectra for 
both the simulated and the experimental system. Several 
singularities can be seen, the most important one is the 
jump singularity on the left edge of the spectrum. If 
fcj — > 0, i.e., if the substrate potential is switched off, 
this singularity shifts to A — and the phonon spectrum 
G(A) takes a finite value at A = 0. This will then lead 
to a divergent integral in eq. (16 1 and, as a result, c(t) 
will logarithmically diverge [25 . This is the well-known 
instability of 2D crystals which we here see to have been 
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FIG. 5: Insets: phonon spectra, as obtained from apply- 
ing eq. (17 1 to the simulated and experimental bandstruc- 
tures X(qj) presented in Fig. [3] Main figures: particle 
auto-correlation functions obtained from the phonon spec- 
tra through a Laplace transformation (solid lines) and by di- 
rect evaluation and averaging of {u na (t + r)u na (t)) following 
eq. (141 (open symbols). A: simulation, B: experiment. 




cured thanks to the stabilizing effect of the substrate and 
the light-induced shift of the jump singularity of G(A) 
away from A = 0. 

The Laplace transforms of the phonon spectra, 
are given as solid lines in the main figures 
51 and are compared to the c(r) obtained from 
processing particle displacements of 6000 (experiment) 
and 2000 (simulation) different configurations containing 
1017 (experiment) and 563 (simulation) particles. For 
the simulation, being based exclusively on the Langevin 
model, we must find perfect agreement between both 
curves, provided our scheme is consistent and the data 
analysis tools are free from flaws. And, indeed, we do 
find consistency of the data. As for the experimental 
system agreement can only be expected as long as the 
Langvin model applies. We observe a reasonable agree- 
ment for the first five seconds (up to r' = 2) and marked 
deviations towards negative values between r' = 2 and 



t' = 10. These deviations are not produced by poor 
statistics, or by the instability of 2D crystals which we 
have supressed by the substrate, but are a real feature 
of the system. Very likely, it is - in a real-space presen- 
tation - the effect of the wave-vector dependence of the 
friction-coefficients. 



IV. DISCUSSION AND CONCLUDING 
REMARKS 



The novel aspect of the present paper is the exper- 
imental determination of the wavelength dependence of 
the friction coefficients in 2D colloidal crystals and, based 
on these data, a comparison of the colloidal dynamics 
to the more simplified Langevin dynamics. Determin- 
ing the tf-dependence of friction coefficients requires to 
measure two quantities independently, the clastic band- 
structure A(tfj) and the phonon decay times T(qj). The 
friction coefficients then follow from multiplying both 
quantities. Such a measurement is possible only with 
video-microscopy data - a technique whose most appeal- 
ing property is that it allows us to directly "see" the 
overdamped dynamics of the colloids, in the form of the 
exponential decay of the phonon auto-correlation func- 
tion in Fig. [2j 

Our main finding is that in the long- wavelength limit, 
the friction coefficient vanishes, i.e., 7(9}) — > for q — > 0. 
This is clearly the case for the longitudinal modes while 
the extrapolated curves for the transversal modes show 
a small offset at q — which however is probably not 
a meaningful extrapolation. The vanishing of friction 
of long wavelength phonons is a very reasonable result. 
Hurd et al. pointed out that in 3D colloidal crystals back 
flow will give an extra damping to longitudinal modes 
such that for q — > the damping is even larger than 
the Stokes friction coefficient. Transverse modes, on the 
other hand, have been shown not to be affected by a back- 
flow damping; the flow induced by the collective motions 
of the spheres now adds constructively, leading to van- 
ishing friction coefficients at long wavelengths. Quite the 
same mechanism can be postulated to work for our 2D 
system: the transverse modes in the limit q — > will 
also induce constructive interference of the flow, an over- 
all flow field occurs going in the same direction than the 
moving colloids and, as a result, the colloidal particle 
show no longer motion relative to the surrounding fluid, 
leading to a vanishing of the friction and thus a faster de- 
cay of the overdamped phonons. For longitudinal modes, 
our 2D system differs from the 3D system in that the 2D 
layer of particles is coupled to a third spatial dimension 
into which the flow is free to move. This third dimension 
prevents the back-flow damping, and a favorable interfer- 
ence of the flow into the third dimension should again be 
the reason for the vanishing of the friction coefficient for 
q — ► 0. However, a more detailed explanation - as well as 
a theoretical reproduction of 7((fj) in Fig. [3]- has to wait 
for a proper hydrodynamic theory describing our system. 
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Such a theory must also account for hydrodynamic effects 
resulting from the walls confining our systems |33j . Our 
results might also be interesting to be reproduced apply- 
ing more advanced simulation techniques [34] taking into 
account many-particle hydrodynamics. 

With these results in mind, it is clear that taking a 
constant friction coefficient is a rather gross assump- 
tion of the overdamped Langevin model. In addition, 
this model ignores that vanishing friction coefficients in 
the long wavelength limit necessarily permits phonons to 
start propagating (though being still damped, but not 
overdamped). Systematic deviations from the exponen- 
tials in Fig. [2] might be seen as a first indication for the 
existence of such modes. However, it is also evident from 
Fig. g] that given these deviations indeed signalize the 
onset of propagating phonons, they are still rather weak 
and far from being able to dominate the overall dynamics 
of the system. The success and failure of the overdamped 
Langevin model is best summarized by Fig. [5] showing a 
particle auto-correlation function c(r) that agrees to the 
Langevin prediction only within the first few seconds. 
Beyond that time, clear devations towards negative val- 
ues in c(r) (anti-correlated behavior) can be observed, 
a feature which according to our previous remarks must 
be due to the combined effect of the wave-length depen- 
dence of the friction coefficient and, possibly, propagat- 
ing modes. Such a time-delayed anticorrelation has also 



been observed in the two-particle experiments and has 
been interpreted in terms of the standard Oseen tensor 
hydrodynamic coupling [18]. Zahn et al. also showed 
that the Oseen term in 2D colloidal suspensions can lead 
to an increase of the self-diffusion [35] . Again: how these 
results connect to our findings and how exactly hydrody- 
namics produces the anticorrelation, can only be clarified 
with an elaborate hydrodynamical theory. It is our hope 
that the present paper can stimulate the interest in such 
theoretical work. 

As well as motivating research on hydrodynamics our 
work might open novel perspectives in studying crystals 
which exhibit non-overdamped particle dynamics. Such 
crystals can be already realized using a dusty plasma [36] 
where the micron-sized dust particles can be observed 
with videomicroscopy similar to colloidal crystals. 
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